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ABSTRACT 


The  AKC  Inventory  Research  Office  (IRO)  has  been  involved 
in  the  development  of  multi-echelon  inventory  models  for  Amy 
applications  for  the  past  several  years.  During  research  for 
a  study  for  the  Joint  Logistics  Review  Board  (JLRB)  by  IRO, 
several  significant  advancements  were  made  in  multi -echelon 
modelling.  Subsequently,  several  more  were  noted  by  TRO  in 
the  literature  and  were  put  to  use.  This  thesis  summarizes 
the  inulti-eche?,on  work  done  at  TRO  with  particular  emphasis 
cn  the  efforts  which  originated  v*  th  and  followed  the  JLRB 
study.  Philosophy  of  snr lysis  is  emphasized  more  than 
mathematical  derivation.  The  techniques  are  exact  w- ere 
possible. 


Introduction 


In  the  past  several  years  the  AMC  Inventory  Research  Office 
i[T.RO)  has  been  involved  in  the  development  and  application  of 
iculti-echelcn  inventory  models.  During  that  time,  several 
reports  were  printed  describing  some  of  the  multi-echelon  work, 
but  until  now  much  of  the  recent  efforts  had  not  been  published 
This  thesis  is  a  summary  and  reference  for  all  of  IRO's  rele¬ 
vant  multi-echelon  work. 

Emphasis  will  be  placed  on  the  analytical  philosophy  be¬ 
hind  the  development  of  the  models  with  the  intention  of 
stimulating  interest  in  the  multi -echelon  area.  Those  models 
documented  elsewhere  ate  given  only  cursory  treatment,  while 
undocumented  models  are  developed  thoroughly. 

Models  of  continuous  review  inventory  systems  will  be  the 
only  ones  considered,  since  periodic  models  have  little  applica 
tion  within  the  Army.  An  appropriate  beginning  is  with  the 
coincident  bu-;  independent  development  of  two  identical  models 
the  AMMIF  mod"!  at  IRO  and  the  METRIC  model  at  RAND.  A 
chronological  history  from  this  point  is  developed  in  this 
report  which  leads  to  a  description  of  all  of  IRO's  current 
multi-echelon  models. 

Sections  2  and  3  cover  the  initial  multi-echelon  work 
which  was  done  in  the  mid  to  late  sixties.  Sections  4  and  5, 
which  are  more  detailed  than  the  others,  cover  some  of  the 
more  recent  work  which  originated  from  IRO  research  during  a 
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study  for  the  Joint  Logistics  Review  Board.  A  model  developed 
at  RAND  in  which  an  error  was  noted  and  corrected  at  IRC'  is 
briefly  covered  in  section  6,  since  the  analysis  itself  is 
a  significant  development.  In  section  7,  a  discussion  of  a 
heuristic  multi-echelon  model  for  the  SAFEGUARD  ABM  is  given. 
2.  The  AMMIP-METRIC  Model 

The  term  "AMMIP-METRIC  model"  will  be  used  to  signify  that 
the  two  are  essentially  identical  models  and  can  be  discussed 
as  one.  The  METRIC  model  has  been  formalized  at  RAND  into 
a  marketable  computer  program  package,  while  1R0  has  kept 
the  associated  c  mputational  techniques  in-house.  Neverthe¬ 
less,  apart  from  computational  techniques,  the  two  models 
have  remained  identical. 

The  best  references  are  [10]  and  [11].  From  a 
theoretical  viewpoint  [11]  is  superior,  and  would  be  more 
valuable  to  one  wanting  to  learn  about  the  model. 

2.1  Basic  Methodology 

A  queuing  theorem  due  to  Palm  is  the  basis  of  AMMIP- 
METRIC  .  Palm  derived  the  distribution  of  the  number  of 
customers  in  an  M/gA>  queue.  His  theorem  states  that  if  T 
is  the  average  service  time,  and  X  the  customer  arrival 
race,  then  the  number  of  customers  in  the  queue  is  Poisson 
distributed  with  parameter  XT,  independent  of  the  form 
of  the  service  time  distribution. 
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A  clear  and  concise  proof  of  Palin's  theorem  is  given  in 
[71  as  a  sidelight  to  another  theorem. 

Often,  inventory  systems  are  studied  by  analogy  with 
queuing  systems.  The  infinite  channel  system  is  frequently 
used.  Demands  are  analogous  to  queue  customers,  and  lead 
times  are  analogous  tc  queue  service  times. 

In  the  multi-echelon  context,  a  lower  echelon  stockage 
location  requisitions  from  the  echelon  above  (queue  customer 
arrives)  and  receives  his  stock  a  lead  time  later  (queue 
service  time).  Now  the  lead  time  may  be  thought  of  as  two 
segments  -  a  normal  time  to  respond  and  a  delay  incurred 
if  no  stock  is  available.  The  lead  time  is  the  sum  of  these 
two.  In  general,  the  delay  due  to  stockout  will  depend  upon 
the  demand  pattern  and  the  stockage  policy  at  tne  above 
echelon.  But,  using  the  queuing  analogy.  Palm's  theorem 
says  that  provided  demand  on  tne  above  echelon  is  Poisson 
distributed,  the  number  of  requisitions  from  the  lower 
echelon  unit  which  have  not  been  filled  (in  the  queuing 
system)  is  Poisscn  and  depends  only  on  the  average  lead 
tim.;,  including  stockout  delays.  Thus,  if  T  is  the  average  of  the 
n-  rrnal  time  to  satisfy  a  requisition  when  there  is  no  stockout  delav 
and  W  is  the  stockout  delay,  the  average  lead  time  is  T  =  Tq  f  E(W). 

At  this  point,  yet  another  queuing  analogy  is  used.  In 
simple  terms,  for  sake  of  discussion,  when  the  lower  echelon 


3 


requisition  arrives  at  the  above  echelon,  it  is  processed 
and  sent  to  the  stock  room  for  shipment.  At  the  stock  room, 
it  is  serviced  and  shipped  without  delay  if  stock  is  available, 
but  waits  at  the  stock  room  as  a  backorder  to  be  serviced 
in  a  FIFO  priority  if  stock  is  not  available.  Here  the  well 
known  L  =  AW  queuing  relation  applies.  Thus,  E(W)  =  expected 
backorders  (L)  divided  by  the  demand  rate  (X). 

Now  suppose  that  all  stockage  locations  follow  S-1,5 
policies  and  that  all  exogenous  demands  are  Poisson  dis¬ 
tributed.  The  S~1,S  policy  has  two  ramifications.  First, 

S-1,S  policies  merely  pass  on  demands  to  the  above  echelons 
with  no  modification  to  their  distribution.  Thus,  S-1,S 
policies  insure  Poisson  demand  at  all  locations  in  the  multi¬ 
echelon  system.  Secondly,  with  S-1,S  policies,  the  net 
stock  at  a  stockage  point  is  S  minus  the  number  on  order. 
Knowledge  of  the  number  of  units  on  order  is,  therefore, 
equivalent  to  knowledge  of  net  stock,  which  is  useful  in 
forming  cost  or  performance  expressions. 

Consider  the  two  echelon  situations  previously  discussed. 
A  stockage  location  in  the  above  echelon  sees  some  average 
response  to  its  requisitions  and  by  Palm’s  theorem  the 
number  on  order  by  the  stockage  location  is  Poisson.  On  the 
average  it  backorders  (net  stock  less  than  zero) 

-CD 

3=2  n*p(n) 
n=I 
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where 

n  *  net  stock  =  on  hand  -  backorders 
p(n)  =  probability  net  stock  equals  n 
=  probability  on  order  equals  S-n 
Then  E(W)  =  B/X  where  X  is  the  demand  rate  on  the  above 
echelon  location.  This  is  then  added  to  the  normal  response 
time  provided  to  the  lower  echelon  and  Palm's  theorem  is  used 
again  to  get  the  distribution  of  net  stock  at  the  lower 
echelon  location. 

Note  that  the  use  of  Palm's  theorem  is  only  an  approxi¬ 
mation.  Delay  at  the  upper  echelon  is  conditional  on  the 
demands  occurring  on  the  lower  echelon.  The  model  does  not 
recognize  this  dependence  and  assumes  that  delay  is  independent 
of  lower  echelon  demands.  Also,  it  does  not  recognize  thcc 
backordered  requisitions  are  not  likely  to  cress  over  in  the 
real  world  as  the  use  of  Palm's  theorem  implies,  i.e., 
lead  times  are  assumed  independent.  Despite  these  faults, 
the  model  has  been  found  to  provide  good  approximations  for 
the  low  demand  items  for  which  it  was  designed. 

2.2  Optimization  Using  A1MIP -METRIC 


There  are  at  least  three  types  of  objectives  for  which 

an  optimization  procedure  is  required.  One  might  wish  to 

minimize  a  total  cost  expression,  or  minimize  the  investment 

required  to  achieve  a  performance  target,  or  achieve  the 

best  perf  unnance  subject  to  an  investment  constraint. 
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Optimization  using  the  AMMIP-METRIC  model  is  made  difficult, 
however,  because,  in  general,  the  objective  functions  are 
not  convex  as  are  those  of  most  single  echelon  inventory 
models.  Often  there  are  small  bumps  in  the  objective 
function  surface  which  prevent  common  optimization  techniques 
that  rely  on  convexity  from  proceeding  to  proper  termination. 
Moreover,  these  bumps  can  arise  with  small  changes  in  para¬ 
meters  with  the  result  that  dramatically  different  alloca¬ 
tions  occur  even  though  the  parameters  differ  by  only  a 
little.  This  was  observed  to  occur  on  a  heuristic  algorithm 
developed  by  IRO.  This  being  a  very  undesirable  property, 

IRO  developed  two  algorithms,  both  of  which  are  considered 
satisfactory. 

The  first  algorithm  is  described  in  [3].  It  was 
designed  to  minimize  total  cost  equal  to  the  sum  of  inventory 
and  backorder  costs  in  a  two  echelon  system.  The  algorithm 
''roduces  exact  optimal  solutions,  but  its  lack  of  applicability 
to  other  than  minimum  cost  objectives  led  IRO  to  develop  a 
flexible  heuristic  algorithm  which  does  not  hang  up  because 
of  non -convexity,  and  which  can  solve  either  of  the  three 
objectives  listed  above. 

The  algorithm  operates  by  adding  one  unit  of  inventory 

at  a  time  to  the  location  where  the  greatest  improvement 

in  total  backorders  occurs.  In  the  terminology  of  search 

procedure,  it  is  a  steepest  ascent  method.  Termination  occurs 
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when  either  tota?  inventory,  stock  availability,  or  total 

backorders  meet  or  exceac  their  targets.  Establishment  of 

an  availability  goal  is  equivalent  to  minimizing  cjsf  (inventory 

plus  backorder  cost}  since  a  necessary  condition  for  minimum  cost 

that  availability  equal  1-C„/C_  where 

n  d 

C  =  bolding  cost  per  unit  per  unit  time 
h 

C_  =  backorder  cost  per  unit  per  unit  time 
B 

In  most  cases,  this  procedure  ..eads  to  the  optimal 
solution.  Where  it  did  not  in  the  cases  examined,  the 
solution  differed  from  optimum  by  no  more  than  one  unit 
at  an/  single  stockage  location.  This  degree  of  error  is 
acceptable. 

3.  Real  Time  Multi-Echelon  Models 

"Real  Time"  as  used  here  denotes  a  stockage  decision 
process  which  relies  on  real  time  system  information  to 
allocate  available  assets  as  opposed  to  allocation  by  a  pre¬ 
determined  decision  rule  such  as  a  reorder  point,  reorder 
quantity  rule  which  does  not  depend  on  any  other  system 
conditions.  There  are  two  models  of  interest  which  are  of 
this  type.  Both  will  be  treated  casually. 

3.1  Real  Time  METRIC 

Reference  [6]  provides  a  complete  description  of  Real 
Time  METRIC.  This  model  provides  a  decision  rule  for  shipping 
depot  stock  to  the  bases  which  is  geared  to  the  occurrence 


of  an  event  in  the  system.  Stated  in  its  simplest  form  the 
rule  affirms  a  shipment  to  base  j  if  bas^  j’s  "need"  is 
greater  than  the  depot's  "reluctance"  to  ship.  "Need"  is 
defined  in  terms  cf  base  backorders,  but  "reluctance"  is 
an  abstract  concept  defined  by  a  parametric  equation,  whose 
parameters  have  been  set  to  yield  the  best  results. 

While  philosophically  appealing.  Real  Time  METRIC  is  somewhat 
lacking  in  rigor. 

3.2  IRQ  Allocation  Model 

In  contrast  to  real  time  METRIC,  this  model,  which  is 
described  in  [2],  operates  only  when  a  stock  imbalance  or 
scarcity  occurs.  A  short  term  horizon  is  defined  and  a 
decision  for  allocation  of  available  assets  is  produced 
which  minimizes  delay  to  customers.  The  model  is  mathe¬ 
matically  sound,  but  is  limited  to  a  three  area  depot,  one 
wholesale  depot  system. 

4.  Delay  Due  to  S:ock-out  at  a  Supplier  With  an  R,Q  Inventory 
Policy 

The  importance  of  delay  due  to  stockout  in  multi-echelon 
models  can  be  appreciatec  rrom  the  discussion  of  the  AMMIP- 
METRIC  model  in  section  2.  This  section  describes  two 
approaches  to  stockout  delay  for  the  more  general  continuous 
review  R,Q  model. 

4.1  Expected  Stockout  Delay 

The  reference  for  this  section  is  [13]  in  which  Simon 
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developed  an  expression  for  average  delay  using  these  assump¬ 
tions  J 

a.  The  demand  process  on  the  supplier  is  Foisscn  with 
parameter  X. 

b.  The  supplier’s  lead  time  is  either  deterministic  or 
exponential.  For  each  lead  time  case,  the  delay  expressions 
were  found  from  the  basic  i*elationship 


E(t)  =  E  E(T|b).?r[B“b] 
b=0 


LA 


where 

r  =  customer  delay  due  to  stockout 
E(r|b)  *  expected  customer  delay  given  b 
backorders  at  his  arrival. 

Pr[B=b]=  probability  backorders  equal  b 

at  the  customers  arrival  (for 
Poisson  customer  arrivals  this 
is  the  same  as  the  probability 
that  backorders  equal  b  at  random 
point  in  time). 

Using  complex  reasoning  involving  order  statistics, 

Simon  was  able  to  find  E(t  |b).  Since  Pr[B=b]  has  been  derived 
in  other  works  (see  [1],  Chapter  4)  he  was  sale  to  find 
E(t). 

His  final  expressions  unfortunately  provided  little 


insight  into  the  delay  process.  However,  for  deterministic 
lead  times,  it  was  shown  in  [8]  that  his  result  simplifies 
to  E[3j/a  which,  of  course,  then  agrees  with  the  intuitively 
appeal  ling  L  =  as’.  However,  Simon  indicates  that  this 
relationship  is  not  correct  for  the  exponential  lead  rimes. 

He  cites  computational  experience  in  which  as  much  as  53% 
difference  was  observed  from  EfBj/A.  This  is  somewhat 
surprising  in  view  of  the  general  applicability  cf  L  =  AU. 

4.2  Probability  Distribution  of  Cu.>  ■‘.omer  Delay 

A  requirement  for  understanding  this  section  is  an  under¬ 
standing  of  chapter  4  of  Hadley  and  Vfnicin  ’ll,  in  which 
the  probabilistic  properties  of  an  R,Q  inventory  oolicy  are 
given.  A  thorough  description  of  the  concents  of  this 
section  is  in  [4]. 

If  the  demand  process  on  the  supplier  is  Poisson  and  if 
his  lead  time  is  deterministic,  then  the  probability  that 
a  customer  arriving  at  time  0  waits  longer  than  r  is 


X  Q 

Hv't)  =  r  U  Pnd(T-X,0)5  R+jj  0  £  T  <:  T  4.2 

W  j=l 


where  d(T~I,0)  =  demand  in  the  interval  [r-T,0 j 

T  =  suppliers  deterministic  lead  time 


This  follows  simply  from  the  fact  that  the  suppliers  assets 
at  r-T,  all  of  which  will  be  available  for  issue  no  later 
the»r  t,  and  the  demand  in  the  interval  [r-T,0)  determines 
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Then 


Var(i-)  =  E(t2)  -  [E(-)J2  ,  (4.5) 

4 . 3  Probability  Distribution  of  Customer  Delay  for  Order 
Size  a  1 

Again  a  complete  reference  for  this  section  is  [4j. 

When  customer  demand  is  always  for  one  unit,  there  is  no 

problem  in  defining  customer  delay.  But  if  a  demand  can  be 

^  1  unit,  then  it  is  possible  that  all,  part,  or  none  of 

the  demand  will  bedelay«ri  due  to  stockout.  As  such,  several 

possible  definitions  of  delay  can  be  made. 

In  order  to  overcome  this  definition  problem,  the 

delay  distribution  is  derived  for  ail  individual  units  cf 

the  demand.  If  a  demand  is  for  U  units,  each  is  identified 

by  an  index  j.  Using  arguments  similar  to  those  in  section 

th 

4.2,  the  probability  that  the  j  unit  waits  less  than  t  is 

»0,j>R*Q,0sTvT 

R-+Q 

G.(r)  =  £  Pr[A'T-T)=a]l’r[d(T-T.O)  a-j  j  (4.6) 
J 

a=max( j,R+l) 

1  <;  j  R-K},  0  t  <  T 


Gj (T)  =  1  tor  all  j  . 

Here  again,  as  with  equation  4.2,  G^(t)  is  strictly  valid  only 
for  a  compound  Poisson  demand  process. 

Derivation  of  the  pdf  by  individual  units  provides 
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flexibility  in  developing  measures  of  delay.  For  example, 
to  get  the  expected  value  of  the  average  wait  of  the  demand 
taken  over  all  units  use 

co  U 

E  (t)  =  2  fiuX  2  E(t.)  . 

url  u  j=l  3 

where 

til 

E{tO  =  expected  delay  of  the  j  unit 
f(u)  =  probability  the  order  size  is  u. 


replenished  in  a  deterministic  lead  time:,  L,  the  probability 
distribution  of  its  net  stock  position  can  be  derived. 
Defining  net  assets  at  time  t,  A(t),  as  on  hano  +  on  order 
-  due  out,  and  noting  that  all  on  orders  at  t-L  will  have 
been  received  into  on  hand  stock  by  rime  L  we  :.uve 

N(tfL)  =  A(t-L)  -  d{t-L,t)  U'.i.! 

where 

N(tju)  =  net  stock  at  time  t  with  lead  time  L 

d(t-L,t)  =  demand  in  the  interval  [t-L,tj. 

Exact  solutions  for  the  steady  state  pdf  of  net  stock  have 
been  found  only  for  compound  Poisson  demand  processes,  but 
successful  approximations  have  been  used.  For  example, 
approximating  the  pdf  of  assets  as  uniform  equal  to  1/Q 
which  is  true  only  when  all  demands  are  of  unit  order  size, 
or  using  a  normal  distribution  for  lead  time  demand. 

Except  for  exponetia ily  distributed  lead  times,  no 
exact  method  exists  for  finding  the  pdf  of  N  for  random 
lead  times.  Had lay  and  Whit in,  however,  suggest  as  an 
approximation 

Pr[K(t)=nj  =  “  Pr'N(t jl.)=n]  g(L)dL  (5.2) 

L=0 

where  g(L)  is  the  pat  of  the  lead  time.  This  will  be  a  good 
approximation, provided  the  chances  of  more  than  one  order 
outstanding  are  negligible.  Equivalent  to  (5.2)  is  the  use 
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of  the  marginal  lead  time  demand  distribution  directly  in 

(5.1). 

If  the  marginal  distribution  of  demand  cannot  be  found 
in  tractable  form,  it  might  be  necessary  to  hypothesize  a 
reasonable  form  for  the  distribution  and  set  its  parameters 
appropriately.  In  this  light  a  useful  result  is  found  in 
Parzen  [10]  where  it  is  shown  that 


Var[X]  =  E[Var(xjY)]  +  Var[E(x|Y)] . 


(5.3) 


Using  lead  time  demand,  d^,  in  place  of  X,  and  L  m  place 
of  Y  in  (5.3)  yields 

Var (d  >  =  E[Var(d  J  L)  j  +  Var[E'd ..II)]  (5.4) 

L  4J 9 

For  compound  Pci a sou  demand  distributions  where  Var  [d.j L]  *  XL  VMR 
and  as  an  approximation  for  others,  (5.4)  car  be  changed  to 


Var[d  ]  =  WR  X  S  E(L)  +  US)'  Var(L) 


(5.5) 


where 


X  =  demand  rate 


S  =  ave.age  order  size 


VMR  =  variance  to  mean  ratio  of  lead  time  demand 
quantity 

Along  with  the  expected  lead  time  denana,  X  S  £(L),  Var(d^) 
can  be  used  to  set  the  parameters  of  the  hypothesized  dis¬ 
tribution  (provided  of  course  it  i-,  a  two  pa-areter  distribution). 
The  military  services  plan  to  determine  their  lead  time 


demand  in  this  manner  using  a  convenient  approximation  to 
the  normal  distribution  which  gives  closed  form  expressions 
for  the  optimal  parameters. 

5.2  A  Two-Echelon,  R,Q  Model 

The  logic  behind  ».nis  model  is  much  the  same  as  the 
AMMIP-METRIC  model,  although  a  few  more  approximations  must 
be  made.  Recall  that  AMMTP-METRIC  developed  measures  for 
the  too  echelon  independent  of  the  lower  echelons,  and 
then  used  these  to  determine  the  effect  of  top  echelon 
stock  on  Lhe  lower  echelon.  In  this  sense,  the  two-echelon 
P.,0  model  is  like  AMMIP-METRIC. 

In  this  case  the  mean  and  variance  of  customer  delay 
due  to  stockout  at  the  top  echelon  are  determined  as  a 
function  of  its  stockage  policy.  Then  these  are  related  to 
the  marginal  lead  time  demand  distribution  to  determine  the 
pdf  of  net  slock  at  the  lower  echelon  locations. 

As  previously  mentioned,  the  use  of  equation  (5.2)  on 
(5.1)  is  identical  to  using 

N(t)  =  A(t-L)  -  dL  (5.6) 

where  d^  has  the  probability  function 

■JO 

Prr d  =X]  =  f  Pr[d(t-L,t)=X]  g(L)dL 
L  L=G 

Consider  a  two  echelon  supply  system  with  several 
stockage  locations  in  the  bottom  echelon  and  only  one  in  the 


top  echelon  from  whom  the  botrom  echelon  points  order. 

As  with  the  AMMIP-METRIC  model,  the  bottom  echelon  lead 
time  is  thought  of  in  terms  of  a  normal  response  plus 
a  stockout  delay  at  the  top  echelon.  In  the  simplest 
case,  the  normal  response  is  deterministic,  with  all 
randomness  coming  from  the  stockout  delay.  Clearly,  if 
g(L)  can  be  found  in  terms  of  stockout  delay,  the  basis 
for  a  two-echelon  model  is  created. 

If  the  lower  echelon  lead  time,  L,  is  equal  to 

C4W 

where 

C  =  deterministic  normal  response 
N  =  delay  due  to  stockout 

then  g(L) ,  L  2  C  is  equal  to  the  stockout  delay  density 
function  h(»)  at  the  point  L-C,  i.e.  g(L)  =  h(L-C). 

However,  even  in  the  simplest  situation  we  were  unable  to 
obtain  a  closed  form  expression  for  the  marginal  distribution 
of  lead  time  demand.  This  was  tried  by  using  a  Poisson 
distribution  to  represent  both  demand  on  the  top  echelon  (it 
cannot  be  if  the  bottom  echelon  locations  order  quantities 
greater  than  1),  and  demand  on  the  bottom  echelon. 

We  decided,  therefore,  to  assume  a  form  for  the  marginal 
distribution  and  set  its  parameters  as  discussed  in  Section 
5.1.  The  negative  binomial  distribution  was  selected  since 
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it  is  gaining  acceptance  within  the  military  services  to 
represent  demand  likelihood.  Of  course,  any  other  distribu¬ 
tion  can  be  used  in  place  of  the  negative  binomial. 

With  L  =  C-W  we  have 
E(L)  =  CHf:(W) 

Var(L)  =  Var(W) 

Using  equation  5.5 

Var[d  ]  =  VMR  AYD  (C+E(W»  +(AYD)2  Var(W) 

Xj 

E[dL]  =  AYD  (C+E(W))  (5.7) 

where  AYD  is  annual  yearly  demand  and  C  and  W  are  expressed  in 
years . 

Demand  on  the  top  echelon  will  depend  on  the  demands 
on  the  lower  echelon  stockage  points  and  their  reorder 
quantities  as  well.  In  section  4  it  was  indicated  that,  at 
best,  measures  of  the  mean  and  variance  of  stockout  delay 
could  be  obtained  exactly  only  for  compound  Poisson  demand 
distributions.  Moreover,  we  have  been  unsuccessful  in  finding 
tractable  expressions  for  anything  but  a  pure  Poisson  demand 
process.  (Equations  4.3  and  4.5)  By  tractable  expression 
is  meant  one  which  can  be  quickly  evaluated  by  computer. 

We  decided,  therefore,  on  the  following  intuitive  approach. 

In  order  to  limit  the  number  of  computations,  the 
reorder  quantities  at  the  bottom  echelon  were  assumed  to 
be  Wilson  Q’s.  This  eliminates  searching  for  the  optimum  Q's 
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Since  the  optimum  will  be  larger  than  the  Wilson,  and  since 
the  tendancy  within  the  services  is  to  keep  order  quantities 
small,  this  is  not  felt  to  be  a  serious  limitation. 

Establishing  the  Q's  establishes  the  demand  pattern 
on  the  top  echelon.  We  will  assume  that  the.  demand  process 
on  the  bottom  echelon  is  Poisson.  Then  demands  on  the  top 
echelon  from  a  particular  lower  echelon  unit  occur  with 
Gamma  distributed  inter  arrival  times. 

That  is  if 

=  demand  rate  on  location  i 

Qi  =  reorder  quantity  of  location  i 

=  time  between  placement  of  orders  on  top 
echelon  by  location  i 

then  is  distributed 

f(u)  =  x. 

r(Q.)  (5.8) 


aiV>Qi"1 


Pelczynski  [9]  has  derived  relationships  for  the  mean  and 
variance  of  the  number  of  order  placements  in  a  random  time 
interval,  n^t),  for  Gumma  distributed  time  between  order 
placements.  In  terms  of  the  parameters  of  (5.8)  then 


1 


The  expression  for  Var[n  (t)]  is  not  as  simple,  but 
much  difficulty  it  can  be  evaluated  o.  a  computer, 
the  limiting  form  as  t  —  >»  is  simply 


without 
However , 
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1 


1  (Q  -1)(Q  +1) 

lim  Var[n  (t) ]  =  — [\.t  +  - - -  ] 

t-x»  L  Q,  1  6 


which  can  be  used  for  some  two  echelon  systems.  In  general, 
the  appropriateness  of  .’is  form  depends  on  the  magnitude 
of  exp(-A.^();  the  smaller  the  better.  In  particular,  if 
the  top  echelon  is  the  wholesale  level,  then  t  would  be  its 
procurement  lead  time.  For  any  but  the  most  inactive  items, 
the  approximation  will  be  good.  Assuming  that  demands 
from  the  lower  echelon  units  are  independent  of  one  another, 
then  the  mean  of  the  quantity  demanded  on  the  top  echelon 
in  a  random  period  t  is 


N  Q  A  N 

E[d(t)]  «  Z  *  2  K  t  (5.9) 

i=l  ^i  i=l  1 

and  the  corresponding  variance  is 

N 

Var[d(t)]=2  Q  Var[n(t)]  (5.10) 

i=l  i  1 

where  N  is  the  number  of  stockage  locations  in  the  lower 
echelon. 

Here  a  critical  assumption  is  used.  While  neat  forms 
for  expectation  and  variance  of  stockout  delay  were  obtained 
for  the  Poisson  only,  we  assume  the  form  of  the  expression 
is  valid  for  any  probability  distribution.  Thus,  wherever 
a  Poisson  probability  function  appears  in  the  expression, 
the  corresponding  function  for  another  distribution  is  used. 
This  we  assume  provides  a  good  approximation. 
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This  is  analogous  to  Hadley  and  Whitin  [1]  replacing 
the  Poisson  by  the  normal  in  the  net  stock  probability 
equation  even  though  the  equations  were  exact  only  for  the 
Poisson.  Moreover,  it  is  reassuring  to  note  that  the  Poisson 
cumulative  distribution  appears  explicity  in  the  expressions, 
and  replacing  the  Poisson  by  another  distribution  does  not 
destroy  the  interpretation  of  the  expression-.. 

Thus,  expected  delay  will  still  be  E[B]/\  as  equation 
(4.3)  was  interpreted.  While  no  intuitive  interpretation 
was  made  for  the  expression  for  variance  of  delay,  any 
interpretation  of  equation  (4.4)  will  not  change  with  a 
substituted  distribution. 

All  the  ideas  having  been  covered,  the  computational 
aspects  of  the  model  will  be  summarized  in  instruction  form. 

1.  Compute  the  for  each  lower  echelon  location  using 
the  Wilson  formula  and  select  reorder  points,  R.. 

2.  Use  equations  (5.9)  and  (5. 1C)  to  determine  mean 
and  variance  of  quantity  demanded  on  the  top  echelon  during 
its  lead  time. 

3.  Assume  demand  on  the  top  echelon  in  the  lead  time 
is  distributed  as  a  negative  binomial  random  variable  and 
compute  its  parameters  using  the  results  of  step  2. 

4.  Assume  a  top  echelon  R  and  Wilson  Q  and  use 
equations  (4.3)  and  (4.4)  with  the  negative  binomial  of 
step  3  to  get  expectation  and  variance  of  stockout  delay. 
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5.  Use  equation  (5.4)  to  obtain  mean  and  variance  of 
lower  echelcn  lead  time  demand  for  each  lower  echelon  unit. 

6.  Use  expression  (5.6)  to  determine  probability 
functions  for  net  stock  at  each  lower  echelon  stockage 
location,  and  also  for  the  top  echelon. 

7.  Form  an  appropriate  objective  function  and  find  the 
optimum  reorder  points. 

6.  An  Exact  Two-Echelon  Model 

In  [14]  a  two-echelon  model  was  developed  and  was  claimed 
to  be  exact.  However,  there  was  an  error  in  the  development, 
Nevertheless,  the  basic  methodological  approach  was  valid. 

A  corrected  methodology  was  developed  in  [5]. 

6.1  Methodology 

The  model  is  developed  from  these  assumptions: 

a.  The  top  echelon,  or  depot,  uses  an  R,Q  continuous 
review  replenisnment  policy  and  a  repair  as  received  repair 
policy. 

b.  The  lower  echelon  locations  use  S-1,3  replenishment 
policies  and  repair  as  received  repair  policy. 

c.  The  demand  process  on  the  lower  echelon  is  Poisson. 

d.  All  repair  and  replenishment  times  are  deterministic. 

e.  There  are  probabilities  that  a  failed  item  can  be  repaired 
locally,  or  if  not  locally  repaired  then  at  depot. 

It  is  beyond  the  scope  of  this  report  to  reproduce  the 


mathematical  expressions,  but  the  ideas  behind  these  are 
interesting  and  will  be  presented. 

When  a  stockage  location  uses  an  S-1,S  inventory 
policy,  knowledge  of  the  number  of  units  on  order  plus 
in  repair  is  equivalent  to  knowledge  of  the  net  stock, 
since  net  stock  plus  on  crder  plus  in  repair  is  always 
equal  to  S.  By  assumptions  c,  d  and  e  the  number  in 
repair  at  time  t  2s  Poisson  distributed.  The  number  on  order 
at  time  t  consists  of  those  orders  which  ha^e  not  had  sufficien 
time  to  be  filled,  plus  those  which  have  had  sufficient  time 
but  are  unfilled  because  of  stockout  delay  at  the  depot. 

If  t^  is  the  order  and  ship  time  for  lower  echelon  location 
j,  then  any  demands  placed  on  the  depot  in  the  interval 
(t-U,t)  cannot  be  sacisfied  by  t  and  are,  therefore,  in  the 
on  order  quantity  at  time  t.  Again,  these  demands  are  Poisson 
distributed  by  assumptions  c,  d,  and  e. 

- 1 - 1 - i - ! - 

t-t  -t.  t-k  -t  t-t.  t 

°  J  o  J  J 

FIGURE  6.1 

Figure  6.1  will  be  helpful  for  the  remainder  of  the  discussion. 

Rq  is  the  depot  repair  time,  and  t  is  the  depot  lead  time, 

where  R  s  t  (a  similar  development  is  required  for  R  >  t  ^ . 
o  o  ^  o  o 

Demands  by  location  j  on  the  depot  prior  to  t-t.  but  later 


than  t 
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7 . 1  The  SAFEGUARD  Supply  and  Maintenance  System 


The  SAFEGUARD  Supply  and  Maintenance  System  will  be  a 
two  echelon  structure,  the  sites  at  the  bottom,  and  a 
depot  on  top*  Removal  and  replacement  of  failed  modules 
(called  ORU's  ter  On-Line  Replaceable  Units)  is  the  basic 
maintenance  concept.  The  failed  ORU's,  most  of  which  are 
reparable,  are  either  repaired  on  site  with  probability,  p, 
or  evacuated  to  depot  on  a  direct  exchange  basis.  Of  those 
returning  to  depot,  there  is  a  probability  q  th*  item  cannot 
be  repaired  at  all. 

The  SAFEGUARD  ICP  will  control  the  depot  operations 
of  maintenance  and  supply.  Depot  stocks  are  replenished  by 
procurement  of  an  amount  Q  when  attrition  reduces  depot 
assets  to  the  reorder  point  R.  Failed  items  returned  to 
depot  enter  the  repair  facilities  immediately  with  no 
batching.  Likewise  repair  at  the  site  is  immediate  also. 

7.2  A  Heuristic  SAFEGUARD  Stockage  Model 

A  most  important  requirement  of  the  SAFEGUARD  ABM  is 
that  a  target  availability  be  achieved.  This  objective  precludes 
the  direct  use  of  a  single  item  inventory  model  in  which 
performance  is  measured  in  terras  of  the  item  alone.  At 
least  two  system  availability  stockage  modf Is  have  been 
developed  to  aid  in  provisioning  ORU’s  to  achieve  the  target. 

But  these  are  single  echelon  models  and  can  only  answer 
the  question  of  how  many  ORU's  are  required  on  site  to 
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achieve  the  target  at  minimum  site  spares  cost.  While  it  is 
possible  to  conceive  of  these  models  being  expanded  to 
handle  the  depot  spares  decision  as  well,  it  could  not  be 
done  without  considerable  effort  and  perhaps  loss  of 
computational  feasibility. 

The  output  from  one  of  these  system  availability 
models  is  a  list  of  spare  ORU  requirements  such  that  the 
system  availability  is  achieved  with  the  least  investment 
in  ORU  spares  at  site.  In  producing  this  list,  of  course, 
the  model  had  to  measure  the  per  dollar  impact  of  a  spare 
on  system  availability.  Spare  ORU’s  affect  system  availability 
through  the  replacement  time  which  is  composed  of  a  normal 
segment  that  does  not  depend  on  spares  plus  a  delay  segment 
that  does  depend  on  spares.  The  output  of  the  availability 
model  can.  therefore,  be  interpreted  as  a  list  of  the  most 
economical  tolerable  delays  due  to  ORU  stockout  on  site. 

Most  important  is  that  as  long  as  these  average  delays  due 
to  stockout  are  achieved,  the  system  will  meet  its  goal. 

Now  average  stockout  delay  is  a  supply  measure  which 
has  been  discussed  throughout  this  report.  In  fact,  all 
the  inventory  models  discussed  are  capable  of  approximating 
average  delay  due  to  stockout  at  the  site  as  a  function  of 
site  and  depot  stockage  policies.  All  of  this  suggests  a 
heuristic  multi-echelon  optimization  procedure  which  uses 
a  suitable  multi-echelon  inventory  model  to  achieve  at  least 
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cost  the  most  economical  tolerable  delays  produced  by  the 
system  availability  model. 

Recalling  the  description  of  the  AMMIP  model  in  section 
2,  Palm's  theorem  was  used  to  approximate  the  number  of  items 
in  the  pipeline  where  demand  for  the  item  was  Poisson 
distributed.  But  at  the  SAFEGUARD  site,  for  a  given  ORU, 
the  number  installed  may  be  small.  Consequently,  failures 
cannot  be  approximated  by  a  Poisson  since  they  are  state 
dependent.  However,  there  is  a  theorem  analogous  to  Palm's 
which  gives  the  probability  distribution  of  the  number  of 
items  in  the  pipeline  when  both  failures  and  lead  times  are 
state  dependent.  A  derivation  appears  in  [12].  If  the 
state  of  the  system  is  denoted  by  m,  demands  are  Poisson 
with  rate  X  ,  and  lead  times  have  an  arbitrary  distribution, 
then  the  state  probabilities  depend  on  only  the  >,  's  and 
the  average  lead  time. 

As  with  the  AMMIP-METRIC  model,  average  lead  time  is 
computed  from  a  normal  lead  time  plus  an  average  delay. 

Then  the  site  pipeline  distribution  is  computed  by  the  above 
theorem.  Thus  if 

T  =  the  average  time  to  return  an  unserviceable 
ORU  to  serviceable  condition  on  site  = 
average  pipeline  time. 

f(x)  =  probability  the  number  in  the  pipeline  is  x 
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then 

fix)  =  C  X1X2...ix.1/TX  (7.1) 

where  C  is  a  normalizing  constant.  Note  that  T  -  p  t  -*- 

(1-p)  (t^-W;  where  t  =  average  s? te  repair  time 

t  ~  average  replenishment  time  from  depot 
s 

W  =  average  delay  due  to  stockout  at  depot. 

To  determine  depot  delay,  we  are  forced  to  return  tc  assuming 
depot  demands  are  Poisson.  A  single  echelon  model  described 
in  [15]  is  used  to  determine  average  depot  backorders  for 
a  given  depot  R  and  Q.  Depot  delay  is  then  computed  as 
average  backorders  divided  by  the  depot  demand  rate.  This 
is  an  exact  expression  and  has  been  derived  ir.  the  same 
manner  as  was  equation  (4.3). 

Using  equation  (7.1)  the  average  number  of  site 
backorders  can  be  round  for  any  site  spares  level,  S. ,  and 
any  depot  policy,  R,Q.  Optimization  ever  and  u  with  a 
Wilson  Q  is  accomplished  by  search. 
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